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We evaluate the performance of modern nuclear energy density functionals 
for predicting inner and outer fission barrier heights and energies of fission 
isomers of even-even actinides. For isomer energies and outer barrier heights, 
we find that the self-consistent theory at the HFB level is capable of providing 
quantitative agreement with empirical data. 
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1. Introduction 

Attaining the goal of a fully microscopic theory of fission has been greatly 
helped by advances in the nuclear density functional theory (DFT) and 
recent developments in the field of high-performance computing. A reliable 
microscopic theory of fission based solely on the properties of the underlying 
effective nuclear interactions would offer unprecedented predictive power, 
thereby giving us more confidence in tackling important questions related, 
e.g., to properties of superheavy elements and astrophysical r-process. 

In the recent years, there have been systematic studies of nuclear fission 
in the DFT framework. In particular, the zero-range Skyrme functional, 1,2 
the finite-range Gogny force 3 ' 4 and the relativistic mean-field approaches 5,6 
have been used to compute barrier and half-live estimates. In parallel, the 
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optimization of coupling constants in energy density functionals 7 has been 
revisited with a particular focus on uncertainty quantification and correla- 
tion analysis. In this context, the UNEDF1 parametrization of the Skyrme 
functional 8 was recently developed with fission applications in mind, opti- 
mized to have the realistic deformation properties. 9 

The goal of this work is to assess the predictive power of the UNEDF1 
functional in fission studies by comparing its predictions to experimental 
data such as barrier heights and isomer excitation energies, and bench- 
marking it against alternative models. We briefly outline the nuclear DFT 
formalism for fission in Sec. 2. We then present and discuss our results in 
Sec. 3, before concluding in Sec. 4. 

2. Theory and Background 

We calculate potential energy surfaces through nuclear DFT with con- 
straints, at the level of the Hartree-Fock-Bogoliubov (HFB) mean field. In 
the particle-hole channel, we use the UNEDFI functional and the SkM* 
functional as a benchmark. 10 In the particle-particle channel, we use the 
density-dependent mixed-pairing interaction. 11 All calculations were per- 
formed with a quasi-particle cutoff energy of E cut = 60MeV. The UN- 
EDFI functional prescribes the pairing strengths; for our calculations with 
SkM*, we determined them such that the calculated pairing gaps coincide 
with the experimental odd-even mass differences in 232 Th. 

We explore configurations that explicitly break axial symmetry, to prop- 
erly account for the inner barrier height, as well as reflection symmetry, to 
minimize the outer barrier height. Calculations are performed with the 
symmetry-unrestricted nuclear DFT solver HFODD, 12 which uses a har- 
monic oscillator (HO) basis. We use the lowest 1140 states from 31 de- 
formed HO shells. For each value of the quadrupole moment Q20, the basis 
is deformed according to «2o = 0.05\/Q20, and the basis frequency is set 
to haj = 6.5 + 0.1 x Q2oe~ ' 02 ® 20 ■ This empirical parametrization gives a 
good compromise between accuracy and speed. 13 

3. Fission Barrier Heights 

While fission barrier heights are not observables in the strict sense, they do 
provide important constraints on theories of fission. Our goal is to evaluate 
and compare the ability of theoretical models to predict quantities such as 
fission barrier heights beyond known data. 
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3.1. Comparison of theoretical models 

In a first step, we benchmark the UNEDF1 functional against alterna- 
tive parametrizations of the energy functional (Skyrme SkM* and Gogny 
D1S 14 ) and the microscopic-macroscopic FRLDM 15 models for the inner 
fission barrier height (Fig. 1), the energy of the superdeformed fission iso- 
mer (Fig. 2), and the outer fission barrier height (Fig. 3). 
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Fig. 1. Inner barrier heights of selected even-even actinide nuclei. Experimental 
data 16,17 (black lines with square markers) are compared to predictions of UNEDF1 
(red lines with circles); SkM* (blue lines with triangles); D1S 14 (magenta lines with 
'x'-es); and FRLDM 15 (green lines with plus signs). 
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For the inner barrier height Ea, we note that the three microscopic 
approaches exhibit the same isotopic trend, with a roughly parabolic Ea (N) 
behavior, and differ essentially by the position of the maximum in the U, 
Pu and Cm isotopes. While FRLDM results are similar to the DFT results 
in the heavy actinides, they show marked differences in the lighter systems, 
especially Th and U. Overall, the excitation energy of the fission isomer 
is much more consistent across all theoretical models used here, with the 
exception of the Cm and Cf isotopes where the FRLDM values are as much 
as twice the DFT value (and the empirical data). Most interestingly, we 
note that all the models considered agree better at predicting Eb than 
Ea'- the RMS of the variances about the mean theoretical prediction for 
each element ranges between 1.32 McV (Cm) to 1.88 MeV (Cf) for E B , 
while it ranges from 1.65 MeV (Cm) to 2.74 MeV (U) for E A . This is 
somewhat counterintuitive, as one would think that the difference between 
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Fig. 2. Similar to Fig. 1 but for bandhead energies of fission isomers. The experimental 
data are from Ref. 18 




Fig. 3. Similar to Fig. 1 but for the outer barrier heights. 



the deformation properties of these approaches would be magnified as the 
elongation increases. 



3.2. Comparison with empirical values 

We now turn to the residuals AE = E cxp — E^ with the empirical val- 
ues taken from Refs. 16,17 When assessing these results one needs to bear 
in mind that the empirical uncertainty on the barrier heights is at least 
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0.3 MeV. The residuals for the inner barrier heights are presented in Fig. 4 
for UNEDF1, SkM*, D1S and FRLDM. Both UNEDF1 and SkM* sys- 
tematically overestimate Ea and exhibit a gentle oscillating trend. The 
FRLDM residuals have a downward overall slope with N, which is espe- 
cially pronounced for the U and Pu chains. The D1S residuals have the 
least scatter for inner barrier heights. The residuals for the fission isomer 
excitation energy are shown in Fig. 5. There are but a few measured data 
points, so we simply note that the self-consistent models exhibit less scat- 
ter than FRLDM. Finally, we show the residuals for outer barrier heights 
in Fig. 6. The UNEDF1 functional demonstrates the least scatter, while 
both D1S and SkM* consistently overestimate the outer barrier heights. In 
addition, the FRLDM residuals exhibit a downward slope with N, which is 
similar to the trend for Ea- 
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Fig. 4. The residuals A.E for the inner barrier heights of selected even-even Th (green 
lines with plus signs), U (red lines with inverted triangcs), Pu (cyan lines with 'x'-es), 
Cm (magenta lines with squares), and Cf (yellow lines with circles) isotopes. 



In this comparison with empirical fission barrier heights and fission iso- 
mer energies, the performance of UNEDF1 is on par with that of the other 
self-consistent interactions and the macroscopic-microscopic FRLDM. For 
the four theoretical models discussed in this work, Table 1 displays the 
RMS deviations from experiment for the even-even isotopes. For inner bar- 
rier heights, the D1S functional yields the best agreement, with an RMS 
of 0.709 MeV, but UNEDF1 agrees nearly as well with an RMS of 1.03 
MeV, i.e., these two models are practically equivalent considering the em- 
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Fig. 5. Similar to Fig. 4 but for bandhcad energies of fission isomers. 
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Fig. 6. Similar to Fig. 4 but for outer barrier heights. The additional blue lines with 
triangles mark Ra isotopes. 



pirical uncertainty on barrier heights. The three self-consistent calculations 
yield a similar level of agreement for the fission isomer energy, while for the 
outer barrier heights UNEDF1 is the closest to the empirical values, with 
an RMS of 0.690 MeV. The FRLDM and D1S predictions for outer barrier 
heights are comparable. 

It is to be noted that the inner barrier heights calculated with all micro- 
scopic models are consistently 0.5-1.0 MeV too high for U and Pu isotopes, 
and are often higher than typical values of one-neutron separation energy 
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Table 1. RMS deviations of the first barrier height 
Ea, the fission isomer En, and the second barrier 
height Eg are given for the selection of even-even nu- 
clei considered. 





UNEDF1 
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D1S 


E A [MeV] 


1.03 


1.52 


1.61 


0.709 


E u [MeV] 


0.357 


0.675 


0.351 


0.339 


E B [MeV] 


0.690 


1.13 


1.39 


1.14 



(about 6.5 MeV). Using such numbers directly in reaction models to com- 
pute the fission spectrum, would lead to, e.g., 240 Pu being non- fissile. We 
should bear in mind, however, the large uncertainty affecting the empirical 
data, the simplifying assumptions used in reaction models, and the sys- 
tematic and statistical errors of DFT calculations. On the DFT side, basis 
truncation errors 13 and beyond mean-field corrections can modify the po- 
tential energy surface in a non trivial way. 2,14 If the zero-point energy is 
taken into account according to the method of Ref., 2 for example, the inner 
fission barrier heights can be lowered by approximately 0.5 MeV. Further- 
more, neutron capture leads to compound nuclei in highly excited states. 
The excitation energy E* has the effect of lowering the fission barriers 
as well 19 so one needs to be careful when using the ground-state barriers 
(E* = 0) in estimates of (n,f) cross sections. 

4. Conclusions 

In this survey of fission barrier properties of even-even actinides, we find 
that the recently developed UNEDF1 energy density functional yields pre- 
dictions that agree well with experimental values and are on a par with, or 
better, than predictions of other self-consistent or macroscopic-microscopic 
models. While the inner barrier heights are systematically overestimated, 
one also needs to bear in mind that empirical values are subject to error 
that is at least 0.3 MeV. As fit-observables are selected for future optimiza- 
tions, it will be valuable to consider if additional data such as one-neutron 
separation energy may help to better constrain fission barrier heights. 
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